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Abstract 

Analytical series expansions for the bending of light in the equatorial plane of a Kerr black hole 
are presented in both the strong and weak deflection regimes. It is critical that these are known in 
analytical form so further analysis can be done for predicting different properties of images formed 
in gravitational lensing. Starting with the exact bending angle in terms of the spin parameter, 
we apply a perturbative scheme for rewriting the bending angle as series expansions in terms of 
the impact parameter of the incident light ray. The asymmetry introduced by the black hole spin 
results in spin-dependent shifts in image positions. We apply our results for the case of a galactic 
supermassive black hole to predict angular shifts of relativistic images from the optic axis. This 
would not be possible without the perturbative expansions in the strong deflection regime, only in 
which relativistic images have a chance of being resolved by future telescopes. 
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I. INTRODUCTION 



Observations of the bending of light near the sun in the early twentieth century and 
Einstein rings in deep space Hubble space telescope images in the last decade or so are 
both examples of tests of general relativity in the weak field limit, or more precisely, the 
weak deflection limit. As the light ray skims by, not far from the horizon, however, a 
richer variety of lensing effects are predicted to occur, awaiting observation by telescopes 
with much higher viewing power than what is available today. Bending angle calculations 



l|, |2|, y, 14, [5] for Schwarzschild and Kerr geometries show that as we approach the depths of 
the gravitational potential, multiple looping of the light ray around the center of attraction is 
possible resulting in what are known as relativistic images (see for example, 6|). In the case 
of a Schwarzschild black hole, as the distance of closest approach ro nears the critical value 
of 3m., where m. = GM/c 2 is its gravitational radius, the bending angle has a logarithmic 
form. When expressed in terms of the invariant quantity b' = 1 — b^/b, where b c = Sy^m. 
this strong deflection limit (SDL) of the bending angle is given by 3] 



-7T + log 



216(7-4 v / 3) 



+ 0[b'] + ... 



(1) 



This result is similar, but not identical, to the Darwin logarithmic term this result takes 
into account an improvement to Darwin's result close to the photon sphere. The detailed 
formulation that makes this modification to the strong deflection limit was recently published 
in fl. 

In this paper, we show that the analysis can be applied to the spinning, or Kerr, black 
hole for the case when the light ray stays on the equatorial plane. In order to determine 
the perturbative series for the bending angle, an explicit expression in terms of just the spin 
and mass of the black hole as a function of the impact parameter is needed as the starting 
point. Such an expression was recently derived in 8|. The bending angle was expressed 
explicitly in terms of the black hole mass m. = GM/c 2 and the black hole spin parameter 
a = J/Mc, where J/M is the angular momentum per unit mass of the black hole. With a 
new definition for b' that includes the non-zero spin of the black hole, we define precisely 
the approach to the critical impact parameter, and thus the strong deflection limit (SDL). 

The weak deflection limit (WDL), on the other hand, in both the Schwarzschild and the 
Kerr case are easily defined in terms of the impact parameter as the limit b — > 00. In this 



limit, the weak deflection bending angle series for the Schwarzschild case is given by 

+ ••• (2) 
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A partial generalization of this result for the Kerr case in terms of the coordinate-dependent 
variable r was obtained by Boyer-Lindquist, Skrotskii, Plebanski [h], [lj], fjj and from a 
numerical treatment by Rauch and Blandford |9l. The Kerr WDL series in terms of the 
impact parameter was recently obtained in a completely different context by Petters |L3J. 

ii this paper, we continue to work in the same perturbative framework that was started 
in ~\ and apply it to the case of deflection of rays confined to the Kerr equatorial plane. 
Bending angles, and therefore the position and magnification of images, depend crucially 
on whether the light ray is traversing in the same or opposite direction to the rotation. For 
this reason, throughout the analysis, we carefully keep track of whether the ray orbits are 
direct or retrograde. The perturbative corrections are obtained in a manner very similar to 
that for the Schwarzschild case presented in 0]. 

The need for analytical results for the bending angle cannot be overemphasized given that 
important lensing variables like image positions, magnifications and time delays depend 
crucially on these results. See for many different approaches towards analytical 

results. Perturbative analysis is extremely useful in order to verify, at least to leading 
orders, the predictions of strong deflection lensing. Specific to the Kerr case, we also need to 
be able to glean the contribution to light deflection arising purely from spin. Furthermore, 
higher order relativistic images can be analyzed only in the strong deflection regime, making 
this an important approach to test general relativity beyond the usual classical tests in the 
weak field regime. 

In Section 2, we start with the bending angle result from 8] and definitions of some of the 
variables to be used in the series expansions. We present the weak and strong deflection series 
expansion terms in Sections 3 and 4 respectively, along with numerical plot comparisons of 
these terms with the exact result. 

As a quick preview, we present here two plots that illustrate the main focus of this paper: 
the 2nd-order strong deflection and the 6th-order weak deflection bending angles are plotted 
alongside the exact bending angle for comparison, for direct and retrograde orbits in Figures 
Q] and [2] showing that our series expansions have excellent accuracy; the SDL series and WDL 
series apply nicely in the two regimes. Note that the two assumptions we have made here 
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FIG. 1: The 2nd-order strong deflection and the 6th-order weak deflection are plotted alongside 
the numerically integrated exact formal Kerr bending angle & (in units of it) for direct orbits. 




FIG. 2: The 2nd-order strong deflection and the 6th-order weak deflection are plotted alongside 
the numerically integrated exact formal Kerr bending angle a (in units of it) for retro orbits. 



4 



are: (1) the light ray stays on the equatorial plane, and (2) the spin parameter stays in the 
range < a < m.. As can be seen in the plots, the agreement between the exact and our 
series expansions is better in the direct than in the retro case. 

In section 5 is presented an application of both series to the simplest case of lensing 
geometry in which the source, lens and the observer are in perfect alignment. We show that 
the usual calculation for the Einstein angle that uses the WDL bending angle is applicable in 
the Kerr case as well. However, since there is no contribution to first order, the Einstein ring 
radius and angle are indistinguishable from the static case. In the strong deflection limit, 
on the other hand, the effect of spin on the angular positions relativistic images is relatively 
much greater. This numerical calculation is analogous to the one for the Schwarzschild case 



found in Virbhadra and Ellis 



3- 



We show that the shift in the image position increases 



with the black hole angular momentum. 

II. FORMAL EXACT BENDING ANGLE FOR THE EQUATORIAL KERR CASE 

Let us now consider a light ray that starts in the asymptotic region and approaches the 
black hole, with ro as the distance of closest approach. It then emerges and reaches an 
observer who is also in an asymptotic region. 

We will use the following convenient notation: 

m. (x a 2 

h— — lu s — — and u = — (3) 
r b s mi 

where b s = sb is the invariant impact parameter. The parameter s will be used to keep 
track of the sign of the impact parameter relative to the black hole spin. : s — +1 for direct 
orbits and s = — 1 for retrograde orbits. Quantities that have s as subscript obey the same 
sign convention. For example, uj s takes on the appropriate sign for direct and retrograde 
orbit. Note that in the limit {u s ,u)q — > 0}, we recover the zero-spin Schwarzschild case, and 
in the limit h — > 0, we have the zero-deflection flat metric limit. With a = J/Mc it is also 
convenient to introduce 

a Jc 

° = m7 = GAP' (4) 
as the "normalized" spin parameter. We will limit ourselves to cases where < a < 1, with 
a = being the Schwarzschild limit and a — 1 being extreme Kerr. Next, we define critical 
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h. 



1 - to. 



parameters analogous to the Schwarzschild case in jj]: 

1 + uj s , 3m. 

and r 



We also define the variable 



h sc V r 



h. 



1 - UJj 
1 + 



(5) 



(6) 



with 



1-31^ 



1 - 



a^o , 3m. 

> 1 . 



r / \l + w s / r 
We have introduced these different quantities for the Kerr case, keeping in mind that 

they should go over to those defined in the Schwarzschild case smoothly when a is set equal 

to zero. So, as shown above, as a — > 0, h sc — > 1 and we recover the definition of h! in In 

both cases, h — > at critical, and /i — > 1 as r approaches infinity. 

source image 




observer 



FIG. 3: Thin lens geometry. 



Critical values for the radius and the impact parameters in Kerr geometry are given by: 

(7) 



r sc = 3m. 



1 - 


a 






b sc 




a 


1 + 






b sc 



and 



(b sc + a) 3 = 27m 2 ,(b sc -a). 



(8) 
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From a lensing perspective, we are interested in impact parameters just beyond the critical 
value (SDL) extending all the way to infinity (WDL). We define the dimensionless quantity 
b' as 

b' = 1 - ^ (9) 
b 

where the insertion of the quantity s guarantees that the b' stays between and 1. Note 
that this definition goes over naturally in the Schwarzschild limit: 

l ~— ' 

As will become apparent in later sections, the "normalized" impact parameter b' is extremely 
convenient and natural for describing the range of values from critical all the way to infinity. 

The lensing geometry is shown in Figure [3] where it is assumed as usual that the thickness 
of the lens plane is much smaller than the distances between source and lens D^g. The 
distance (from the observer) to the lens and source are denoted by Di and D$ respectively 
The angular position of the source, (3 and the image, 9 are also shown in the figure. 

The deflection of the light ray from its original path at is given by the following 
expression s| . 



a = -7T + -^— X [^-\Q + [n(n + , k) - IT(n + , V, k)] + fi_ [n(n_, k) - IT(n_, k)] \, (10) 
1 — cu s y Q y J 

where H(n±, k) and H(n±, if), k) are the complete and the incomplete elliptic integrals of the 
third kind respectively. The argument k 2 is defined through the elliptic integral as usual in 
the range < k 2 < 1. Note that in some references the variable is referred to as k 2 and 
in others simply as k. The order in which the arguments appear in U(n,ip,k) also varies 
between different references and in Mathematica. 

Remark: In Mathematica, the built-in mathematical function for the incomplete elliptic 
integral of the third kind EllipticPifn, <ft, m] is defined by 

/ [1 - n sin 2 9] 1 [1 - m sin 2 9] ~ 1/2 d9 
Jo 

and the complete elliptic integral of the third kind is EllipticPi[n, m] = EllipticPifn, it/2, m]. 
In the limiting case when a — > 0, we have f2 + = 1, = and n+ = we have the 
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Schwarzschild result 



-7r + 4 4 /-^[n(0,A;) -n(^,o,fc)] 
-„ + 4J'JL[K(k)-F(il;,k)], 



(11) 



where K(k) and F(ip, k) are the complete and incomplete integrals of the first kind respec- 
tively. In addition, in the limit when m. = (i.e., h — > 0) we recover zero deflection as 
expected. 

The different variables are defined as follows: 



ro 
Q 



k 2 = 



ip = arcsin 



h x 



2h\ / 6h\ 6h 1 



1 - 


2h 

h sc ]j 


K- 


2/i\ 
/W 




1 - 


6h 

h sc \j 


'(- 


2/A 
hsc) 


+ h S c) 



(12) 



(13) 



(14) 



n± = 



±(l±y/T^)(l-u s )^u /2 



LJ 



u h s 



6h 

her 



2h 

her 




^ C V V \ h sc ) u h sc ° l 



(15) 



(16) 



We note here that the quantity b' appears in this expression via r , h, h sc , and u> s , while 
u>o = a 2 /ml is independent of the impact parameter. Any quantity that has an "s" in the 
subscript takes on a negative sign when on the retro side. The angle a itself stays positive, 
i.e., the light ray is still deflected towards the axis of rotation albeit to a lesser extent on the 
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retro side [8|. In the next two sections, we outline the series expansions for the weak and 
strong deflection limits. 



III. EXPANSION OF THE BENDING ANGLE BEYOND CRITICAL 

As we did in [7|, we show that the SDL bending angle can be expressed as an "affine 
perturbation" series in V . The quantity b' is now defined appropriately for the equatorial 
Kerr case. The choice of variables turns out to be very important for carrying out the entire 
calculation successfully. Moreover, the series terms in the equatorial Kerr case become much 
lengthier because of the contribution from the non-zero spin. 

We define an affine perturbation series about a function g(x) as 

f(x) = (A + --- + A p x p + ---)g(x) 

+ (B + --- + B g x q + ■••), 

where A, and Bi are constants with p and q positive rational numbers. Analogous to the 
Schwarzschild case, the the bending angle has an invariant affine perturbation series of the 
form 

a(b') = (a + a 1 (6') + a 2 (6') 2 + ---) log 

(po + Pl (b')+p 2 (b'Y + •••), (17) 

where Ao, o~i and p\ are not just numerical constants, but also depend on the spin parameter a. 
Note that (ITT)) is not a Taylor series expansion because of the appearance of the logarithmic 
term. 

The formal expression ( fTOl) for the bending angle is just that; it is difficult to extract 
information about position and magnification of images as seen from the observer's vantage 
point. The difference between the direct and retro orbits, for example, is simply not at all 
obvious if we just look at the formal exact expression. The effect of frame-dragging on the 
light ray is buried deep in the details of the expression. In what follows, we have sorted 
through these details being careful at each step to make sure that the Schwarzschild results 
are recovered whenever the spin parameter is turned off. 

We begin by rewriting the expression in terms of the variables h, h sc , u s , and u . The 
two regimes that we are referring to as SDL and WDL involve different series expansions of 
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elliptic integral of the third kind and these have to carried out using the built-in EllipticPi 
functions and their properties in Mathematica. The WDL is easily expressed as the limit as 
b — > 0, while a precise definition of the SDL is a bit more involved. Following our scheme 

Q 

in [7|, we introduce the variable b' = 1 — sb sc /b and consider the limit as b' — ► 0. In this 
limit, we are approaching the critical impact parameter. Figure H] is a schematic to show 
the limiting values of the different variables as the impact parameter goes from critical to 
infinity. 
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P = Q = 3m. 







FIG. 4: A schematic to show the strong and weak deflection limits in terms of the relevant variables. 



On the left hand side, we have all parameters and the corresponding limiting values as we 
approach critical radius. First, r is defined in a similar way as we did in the Schwarzschild 
case |7|. The argument k 2 is defined through the elliptic integral as usual in the range 
< k 2 < 1. Note that in some references the variable is referred to as k 2 and in others, 
simply as k. The order in which the arguments appear in H(n, ip,k) also varies between 
different references and in Mathematica. The variables h! and V are both equivalent ways of 
describing the range of possible impact parameters. We will eliminate h' in the end to express 
all quantities in terms of b' in the end. The impact parameter b itself can be positive or 
negative and the sign is carried by s. Note that the impact parameter explicitly approaches 
critical on the left and infinity on the right. Some of the other intermediate variables are 
also given here for completeness. Once again, our goal is to express the bending angle in 
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terms of b f ; the only other quantities that will remain in our final results are m., a and s. 
We have assumed that a stays in the range < a < 1. In addition to the variables shown 
in Figure HI the argument n± also stays in the range {0 — > 1}. 



IV. WEAK DEFLECTION LIMIT 



The weak deflection limit is simply defined by the limit as b — > oo. Series expansions 
for the elliptic integrals are readily available in Mathematica. Before we present our result, 
we give here the existing correction to the Einstein bending angle that other authors have 
derived before. For example, the following result appeared in 1967 in one of the earliest 
papers on the Kerr metric by Boyer and Lindquist 10]: 

& = A^(l + ±) (18) 

The r is the Boyer-Lindquist "radial" coordinate or "the distance of closest approach". 
This result is for direct orbits only and it clearly shows that there is no spin-dependent 
correction to the Einstein angle to first order in 1/6. In the above result and in a number of 
other references, the direct orbit is one in which a is taken to be negative; this is a matter 
of convention. Similar results (some with the opposite sign convention) were also found by 
Skrotskii(ll| and Plebanski[12| and others. Extensive numerical analysis can be found in 
Rauch and Blandford js]. 

As in the Schwarzschild case, we have chosen to keep the expression in terms of the in- 
variant quantity, b instead of tq. After a straightforward calculation and some simplification 
we obtain the following bending angle WDL series for the equatorial Kerr: 
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64 16 

The spin parameter a and the sign s appear explicitly in our series expansion. We will first 
note that the leading behavior of our result is in agreement with (ITS]) after correcting for 
the sign; specifically, we do not see a correction to first order in 1/6. As mentioned in the 
introduction, the leading terms of our series are in exact agreement with those obtained by 
Petters in a completely different context fl3 |. Note that some of the terms that depend 
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on the spin parameter are also sign sensitive. Although it may never be needed, the weak 
deflection series can be extended to include terms of much higher orders, for different values 
of a and for both direct and retro orbits. As a first check of our result, we set a = in 
equation f[T9l 



and recover the Schwarzschild series 

In Figures [TlfTUl plot comparisons of the perturbative and exact Kerr bending angle are 
presented for a = 0.5 and a = 0.99 for both direct and retro orbits. Corrections up to first, 
second and third order in m./6 are plotted with the exact bending angle in each case. The 
plots clearly show that the weak deflection approximation gets closer to the exact angle 
as we include higher order terms. We have shown direct and retro plots to illustrate the 
effectiveness of the weak deflection series for both s = +1 and s = — 1. We have chosen two 
representative values of a for the plots; in fact, the series works nicely for all spins < a < 1. 

The correction to the Schwarzschild weak deflection bending angle coming from spin 
appears in second order (1/& 2 ) only. This means that the spin contribution would be nearly 
impossible to detect in the weak deflection regime. The Einstein ring, for example, when 
the source, lens and observer are perfectly aligned would be essentially the same for the 
Schwarzschild and Kerr case, making it extremely difficult to detect the spin contribution. 
As we will see shortly, the formation of relativistic images in the strong deflection limit is 
much more sensitive to the spin and for this reason would be a far more important tool for 
testing general relativity in the strong field limit and for studying the black hole spin itself. 




V. STRONG DEFLECTION LIMIT 



As we approach critical impact parameter, we once again see, as in the Schwarzschild 
case {2 1, that the leading behavior is logarithmic. In addition to this virulent term, the 
dependence on spin makes the series expansion quite complicated. We have nevertheless 



been able to extract the series terms up to second order in b' as defined in equation (ffj 
The built-in series expansions in Mathematica have been used for this purpose [20!] . Care 
must be taken with defining and providing the correct arguments for the Mathematica 
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functions and appropriate limits for the associated series. We give here only the zeroth 
order term for the bending angle: 



"SDL = 



where 



3W /T 



2^/T^u~ [3 - 2h sc {l - u s )\ log ^ 2(2 ^ — j + 



17_VL + U+V+ 



V 1 - [9 - h sc (6 - /i sc u;o)] (1 - u s ) 



U± = h sc [±u t 2(1 - u s ) (1 ± VT^u~ )] ± 3 [1 ± v / T^ - 2cu s ] 



h sc uo 



6 ± 6\/l — cl> + /i sc u;o 



log 



h sc uo 



6 ± 6\/l — wo + h sc u 



3h sc uo 



6 ± 6\/l — u + /i sc a;o 



h sc Uo 



6 ± 6i/l — a; + h sc u 



3h sc UQ 



6 ± 6y/l — u + /i sc cj / . 



We have written the bending angle in terms of h! because the full version explicitly in 
terms of b' is too complicated to show here. It can be shown, almost just by inspection 
of the above zeroth order term, that setting a equal to zero (which means u s , u — > and 
h sc 1) yields the Schwarzschild version [7| immediately: 

a = -7T + 2 log [ 12{2 ~^ + 0[h'\ + ... (21) 

In Figures [mrHl are shown plots of the strong deflection approximation and the exact 
bending angle. Due to the presence of both positive and negative terms in the series, the 
perturbative expression oscillates about the exact value and gets closer to the exact as we 
include more terms in the expansion. 

VI. EINSTEIN RING AND PARA IMAGES 

The weak and strong series expansions are applicable to either ends of the entire range 
of b' . When both the WDL and SDL series are taken together, they cover the entire range 
with accuracy. Figures [15] and [16] depict this for the case when a = 0.5. The value of b' 
where the discrepancy plots criss-cross varies, depending on the spin value. In Figures [1] and 
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[21 we showed plots of the 2nd-order strong series and the 6th-order weak series along with 
the exact for spin value of a = 0.5. From all these plots it is clearly evident that our series 
expansion results are excellent approximations to the exact bending angle. The is essentially 
the same for all spin values between < a < 1. 

We present here an application of the strong deflection limit of the bending angle used 
along with the lens equation for the lens geometry shown in Figure [3] with the usual as- 
sumption that the thickness of the lens plane is negligible compared to the distances be- 
tween the source, the lens and the observer. For the calculations shown here we will use 
the same model as in [3], where the mass of black hole is taken to be 2.8 x 1O 6 M . With 
Dls/Ds = 1/2 and Di = 8.5 kpc, the angular radius of the Einstein ring can be calculated 
using 6 E = (4m.D LS /D L D s ) 1/2 . The Einstein ring is formed by a cone's worth of rays 
emanating from the source and symmetrically deflected by the black hole. In the Kerr case, 
we expect that the projection of these rays on the sky would still be a closed curve, except 
it won't have circular symmetry; depending on how the magnification varies, these could 
appear as arcs on the sky. Furthermore, our results apply only to light rays that stay on 
the equatorial plane, and hence the plane of incidence. Rays coming in at an angle to the 
equatorial plane, however, will not stay on the plane of incidence: the black hole spin will 
peel these away from the plane of incidence (see [181] and 17J]). Since we are working with 
just the equatorial Kerr case, our results pertain to only two of these rays, one on each side. 
If these were the only two rays emanating from the source, we would see two point images 
on either side of the black hole. These are basically the antipodal points of the "ring" image 
that would be formed. We refer to these image positions as "para images" — "stardogs" 
or "parastars" are terms that would be appropriate as well. Angular positions for the para 
images are listed for the case when the source, lens and observer are perfectly aligned (i.e., 
(3 = 0). Tabled is analogous to Table III in 19|. (The small discrepancies in the a = data 
compared to those obtained by [3] is due to rounding off in the physical constant and in 
conversion factors.) We use the weak deflection limit expansion for calculating 9e for the 
stardog positions on the sky. The angular separation between the para, or stardog positions 
in the Kerr case is analogous to diameter of the Einstein ring in Schwarzschild geometry. 

In the strong deflection limit, i.e., when the impact parameter is close to critical, we see 
multiple loops of the light ray as can be seen from the plot of the exact bending angle. 
Rays that loop around multiple times result in multiple relativistic para images on either 
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Relativistic 




FIG. 5: Para positions or antipodal points on either side of the Kerr black hole where relativistic 
image I would appear on the sky, are shown as a function of increasing a. The data used for this 
hgure can be found in Table U 




FIG. 6: Para positions or antipodal points on either side of the Kerr black hole where relativistic 
image I would appear on the sky, are shown as a function of increasing a. The data used for this 
hgure can be found in Table 21 
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2ir + 43.3477 /ias 


1.93710 


1.01452 




0.9 


-22.2255 nas 


2ir + 44.4510 /ias 


1.98149 


1.01348 




0.99 


-22.7161 lias 


2ir + 45.4321 /ias 


2.02077 


1.01264 


Relativistic Einstein ring II 





16.8996 /ias 


4tt + 33.7992 /ias 


1.50188 


1.00125 


Relativistic Para image II: 


0.000001 


16.8996 nas 


4jt + 33.7992/ias 


1.50187 


1.00125 


(direct) 


0.1 


16.2393 /ias 


4jt + 32.4785/ias 


1.44337 


1.00158 




0.2 


15.5568 fias 


4ir + 31.1135 /ias 


1.38239 


1.00203 




0.3 


14.8477 fias 


4tt + 29.6954 /ias 


1.31850 


1.00265 




0.4 


14.1060 /ias 


4jt + 28.2120 /ias 


1.25112 


1.00356 




0.5 


13.3230 /ias 


4jr + 26.6460 /ias 


1.17944 


1.00494 




0.6 


12.4853 /ias 


4ir + 24.9707 /ias 


1.10226 


1.00713 




0.7 


11.5708 /ias 


4ir + 23.1417 /ias 


1.01765 


1.01091 




0.8 


10.5376 /ias 


4jt + 21.0752 /ias 


0.92208 


1.01826 




0.9 


9.2856 /ias 


4jt + 18.5713 /ias 


0.80701 


1.03606 




0.99 


7.7345 /ias 


4ir + 15.4690 /ias 


0.67348 


1.15357 


Relativistic Para image II: 


0.000001 


-16.8996 /ias 


4jt + 35.7992/ias 


1.50188 


1.00125 


(retro) 


0.1 


-17.5411 /ias 


4ir + 35.0822 /ias 


1.55824 


1.00100 




0.2 


-18.1664 /ias 


4jt + 36.3328 /ias 


1.61272 


1.00082 




0.3 


-18.7774 /ias 


4jt + 37.5549 /ias 


1.66554 


1.00067 




0.4 


-19.3760 /ias 


4ir + 38.7519 /ias 


1.71688 


1.00056 




0.5 


-19.9633 /ias 


4jt + 39.9265 /ias 


1.76687 


1.00047 




0.6 


-20.5405 /ias 


4ir + 41.0811 /ias 


1.81564 


1.00039 




0.7 


-21.1087 /ias 


4ir + 42.2174 /ias 


1.86329 


1.00033 




0.8 


-21.6686 /ias 


4jt + 43.3372 /ias 


1.90992 


1.00028 




0.9 


-22.2209 /ias 


4jt + 44.4418 /ias 


1.95561 


1.00024 




0.99 


-22.7120 /ias 


4jt + 45.4240 /ias 


1.99598 


1.00021 



TABLE I: Para positions of Einstein and relativistic Einstein rings and the analogous images in Kerr 

geometry. Images on the retro side are indicated with negative values for the 9 in microarcseconds 
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(fias ). 



side just as in the Schwarzschild case. The image positions we obtain for direct orbits seem 
reasonable when compared with the higher order images studied by others for the non-zero 
spins for direct orbits off the equatorial plane (see Table I of 18] and [l^] and references 
within.) We believe that image predictions for retro orbits for all spin values < a < 1 have 
not appeared in literature before. Image positions for direct versus retro orbits illustrates 
the key difference in the effect the black hole spin has on rays that are traversing upstream 
versus downstream: the para images are pushed outward on the retro side and inward on 
the prograde side. The data used to generate Figures [5] and [6] are shown in Table fl] Positive 
and negative values of the angular position 9 of the image indicate respectively whether the 
image is on the right or left side. The last column in Table [I] gives the ratio r /r crit , rather 
than r /2m. showing clearly that the rays are just outside critical radii. We would like to 
reiterate here the importance of analytical expansion for the bending angle, without which 
calculation of these image positions would not be possible. 



VII. CONCLUSIONS 



Analytic solutions in the form of series expansions in the strong and weak deflection 
regimes for the bending angle of light rays on the equatorial plane of the Kerr black hole 
were presented in this paper. Higher order terms in strong deflection limit were not shown 
because they are too complicated. The technique to generate these, however, have been 
presented in detail. We have also shown that both expansions have a high level of accuracy 
compared to the exact bending angle, and so would be excellent predictors of image positions. 
To illustrate this, both expansions were applied to the simplest gravitational lensing situation 
when the source, lens and observer are perfectly aligned on the optic axis. We showed that 
the effect of the black hole spin distorts and shifts the Einstein ring inward on the direct side 
and outward on the retrograde side. For large spins, the shift in the image positions is much 
higher than the angular separation of successive relativistic images. Of course, whether 
this shift can be observed with precision will depend on the resolution power of future 
telescopes. The series expansions will be applied next to the case when the source is not on 
the optic axis and also calculations pertaining to magnification of the para images. In order 
to obtain the full geometric structure of images resulting from strong gravitational lensing, 
we plan to pursue a similar analysis of the full Kerr bending angle beyond the equatorial 
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plane. Our goal is to study the dependence of the bending angle on the inclination and the 
resulting variations in the image positions and magnifications in a perturbative framework 
that without any assumptions about the inclination angle or the black hole spin. 
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FIG. 7: Plot comparison of perturbative and exact angle for a = 0.5 direct orbit. 



19 



1= 

•v. 







I I 1 I 


1 1 1 1 1 1 1 1 1 1 

Kerr exact 


- 




WDL order 1 






WDL order 2 




- 


WDL order 3 


- 




(retro, §=0.5) 




1 


i , , T'T*": : — : — i — r— 





0.2 0.4 0.6 0.8 1 

b' 



FIG. 8: Plot comparison of perturbative and exact angle for a = 0.5 retro orbit. 
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FIG. 9: Plot comparison of perturbative and exact angle for a = 0.99 direct orbit. 
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FIG. 10: Plot comparison of perturbative and exact angle for a = 0.99 retro orbit. 
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FIG. 11: Plot comparison of perturbative and exact angle for a = 0.5 direct orbit. 
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FIG. 13: Plot comparison of perturbative and exact angle for a = 0.99 direct orbit. 
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FIG. 14: Plot comparison of perturbative and exact angle for a = 0.99 retro orbit. 
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FIG. 15: Percentage discrepancy for 2nd-order strong deflection and the 6th-order weak deflection 
for direct orbit with a = 0.5. The horizontal axis corresponds to the exact value. 
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FIG. 16: Percentage discrepancy for 2nd-order strong deflection and the 6th-order weak deflection 
for retro orbit with a = 0.5. The horizontal axis corresponds to the exact value. 
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